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Abstract 

Extremal spacings between unimodular eigenvalues of random unitary matrices of 
size N pertaining to circular ensembles are investigated. Probability distributions for 
the minimal spacing for various ensembles are derived for N = 4. We argue that for 
large matrices the average minimal spacing s m i n of a random unitary matrix behaves as 
jy- VC-+0) f or p equal to 0, 1 and 2 for circular Poisson, orthogonal and unitary ensembles, 
respectively. For these ensembles also asymptotic probability distributions of minimal 
spacing s m j n are numerically studied and the statistics of the largest spacing s max are 
investigated. 

1 Introduction 

Random unitary matrices are useful to describe spectra of periodic quantum systems the 
classical analogues of which are chaotic [TJ [2] . The choice of a specific ensemble of matrices is 
dictated by symmetry properties of the investigated physical system. If the system possess no 
time-reversal symmetry the circular unitary ensemble (CUE) of matrices distributed according 
tho the Haar measure of the unitary group is appropriate [3] . For systems with a generalized 
time reversal symmetry the circular orthogonal ensemble (COE) describes properly statistical 
properties of spectra if we neglect additional subtleties caused by specific rotational symmetry 
features of systems with half-integer spin which are of no concern for investigations reported in 
this paper. In the case of classically regular dynamics the spectrum of the evolution operator 
displays level clustering characteristic to the circular Poissonian ensemble (CPE) of diagonal 
random unitary matrices. To describe intermediate statistics one uses interpolating ensembles 
of unitary matrices [4j |6] or composed ensembles of unitary matrices [7] . In the case of 
emerging chaos, in which the chaotic layer covers only a fraction of the phase space of the 
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classical system one may apply the distribution of Berry and Robnik, originally used for 
autonomous systems [8]. 

To characterize statistical properties of spectra of random matrices one often uses the 
distribution of nearest neighbour spacings s, also called level spacing distribution, with den- 
sity as usual denoted by P [31 [3] • The random variable s is the distance between adjacent 
cigenphases normalized to make the mean spacing equal to unity. 

In this work we investigate the distribution of yet another random variable - the minimal 
spacing between two eigenphases s m - m . As the above mentioned distribution of spacings also 
s m in encodes interesting information about the spectrum. Indeed, observe, e.g., that for any 
unitary matrix U the size of its minimal spacing s m i n , provides an information, to which 
extent the matrix U is close to be degenerated. For completeness we are also going to study 
the size of the largest spacing s max defined similarly. 

The distribution of extremal eigenvalues of random Hermitian matrices was analyzed sev- 
eral years ago [TU]. Our current approach is in some sense analogous, as we investigate the 
extremal gaps between eigenvalues of a unitary matrix distributed along the unit circle. After 
a part of our project was completed we learned about a relevant paper of Arous and Bour- 
gade in which the distribution of extremal spacings was studied for random matrices of 
circular unitary ensemble. 

The paper is organized as follows. For exemplary cases of random matrices of size N = 4 
we derive in Section [5] exact forms of minimal spacing distributions. The chosen dimension 
allows exact calculations, which become prohibitively complicated for larger matrices and, on 
other hand, is the minimal dimension in which results for CUE and CPE can be compared 
with those for the ensemble consisting of tensor products of two CUE random matrices of size 
N = 2. Such an ensemble corresponds to a generic local dynamics in a two-qubit system [12] , 

The case of large matrices is studied in Section[3j In particular, we argue that for a random 
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unitary matrix of size N the size of the minimal gap scales as s m - m w N 1 +t> where p = 0, 1 
and 2 for the Poissonian, orthogonal and unitary circular ensemble, respectively. Using known 
formulae for the level spacing distributions we show how to estimate the mean maximal 
spacings for the three ensembles. Finally we formulate plausible conjectures concerning the 
asymptotic (in N) distributions of the minimal spacing s m j n for all three ensembles and 
compare them with numerical results. 

We use the following notation. For a single unitary or orthogonal matrix A of size N we 
shall consider its spectrum {exp(it^j)}^ =1 , where (<pi, . . . , <Pn) represents the vector of the 
eigenphases ordered nondecreasingly, < epi < . . . < ip n < 2tt. We order nondecreasingly the 
spacings ip2 — (pi, ■ ■ ■ , <Pn ~ fN-i, 2tt + fx — <Pn between neighboring eigenphases and denote 
obtained sequence by 

Smin ■= s\< ... < s N =: s max . (1.1) 

The standard level spacing distribution P is given by the average Sm=i where P m is 
the density of the variable s m . 

2 Case study: minimal spacings for TV = 4 

Our goal is to derive the exact probability distribution of the minimal spacing P m i n for ex- 
emplary ensembles of random unitary matrices of size N = 4. Apart of the Poissonian and 
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the unitary ensemble we analyze also the tensor product of two independent random matri- 
ces of size N = 2. This ensemble, denoted by CUE 2 ® CUE2, describes dynamics of two 
independent quantum sub-systems |12j . 

To derive the desired distribution we calculate the tail distribution T(t) = P (s m i n > t) 
and take the derivative of T. We have 



T{t) = P(s min >t) = F(ip 2 - <fi,<P3 - ^2,^4 ~ <<53,27r + ip 1 - ip 4 > nt/2) 

P md (<pi,<P2,(f>3,<PA), 

{<p 2 —<pi,<p3—<p2,<P4—<P3,2Tr+<pi—<p4>wt/2} 



where P ord is the joint probability distribution of ordered eigenphases, which can be obtained 
from the joint probability distribution for a given ensemble. After changing variables, ipi = 
fii "02 = ^2 — yijV's = ~ V2 and ipi = ^4 — V?3i the integration domain splits into two 
tetrahedrons. Simple but long calculations yield in each case the tail distribution function T, 
which leads to the probability density P m ; n = — 4^T. 



(a) CUE 2 ®CXJE; 
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-Pmin(Smin) = ^~ (^(l - S min ) (4 - COs(7TS min /2)) - 3 Sm(7TS min /2) 



+ 8sin(7rs min ) - 3sin(37rs min /2) 



(2.2) 



(b) CUEi 



(c) CPEi 



,(smi„) ==-5 sm z (ns min /4) 666 + 720^(1 - s mm r 

727T Z \ 

+ 36(11 + 16tt 2 (1 - s min ) 2 ) cos(7rs min /2) 

+ 18(8tt 2 (1 - s min ) 2 - 13) cos(7rs min ) - 100 cos(37rs min /2) 

— 608 cos(27rs m i n ) — 380 cos(57rs mm /2) + 234 cos(37rs mm ) (2.3) 

+ 74cos(77rs m i n /2) — 58cos(47rs m in) + 10cos(97rs m i n /2) 

+ 247r(l — s m ; n ) ^60 sin(7rs m ; n /2) + 63 sin(vrs m i n ) + 22sin(37rs min/2) 
+ 2 sin(27rs min ) - 4sin(57rs min /2))j 

^min(Smin) — ^(1 Smin) . (2-4) 



The graphs of these functions are presented in Figure 12.11 The behavior around zero of 
the densities encodes information about level repulsion and level clustering. The variable 
s m in is the smallest distance between two neighboring eigenphases. Therefore, the fact that 
its density is separated from zero, say P m in > 1 around zero, means that for small e > 
the probability that some two phases are in the distance closer than e equals P (s m - m < e) = 
Jo Pmin(s)ds > €. In the cases of CPE\ and CUE4 this consistency with the level clustering 
and repulsion observed in the distribution of spacings P is clearly visible. The figure shows 
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Figure 2.1: Probability densities P m in for the random matrices of size N = 4 pertaining to 
(o) CPE 4 , (□) CUE 2 ®CUE2, (V) CUE4. The symbols denote numerical results obtained for 
2 14 independent matrices, while the curves represent distributions (|2.4[) . (|2.3[) and (|2.2j) of the 
variable s m i n , respectively. 



that for CUE2 ® CUE2 ensemble the cigcnphases of the tensor product tend to accumulate 
in a spectacular contrast to the case a single random unitary matrix form CUE [12] . 

Numerical results show that for large N and m of order of N/2 the distributions of the 
variables s m are close to the distribution of the level spacing s. However, for any N the 
distributions of the smallest spacing s m j n and of the largest spacing s max differ considerably. 
We shall then analyze these distributions of extremal spacings, which can be used as auxiliary 
statistical tools to characterize ensembles of random matrices. 



3 Extremal statistics for large matrices 

In this section we analyze extremal gaps in spectra of circular ensembles of random matrices 
of a large size, N S> 1. As usual, we parametrize canonical ensembles by the level repulsion 
parameter j3, equal to 0, 1 and 2 for Poissonian, orthogonal and unitary ensembles respectively. 
Then, relevant quantities are dentoted with the index j3, e.g. Pp, ft = 0, 1, 2 is the density of 
the level spacing of the corresponding ensemble. 

3.1 Mean minimal spacing 

To get an estimation of the behavior of the mean minimal spacing with the matrix size 
let us assume that spacings Sj, j = 1, . . . ,N are independent random variables. For small 
spacing one has Pp(s) ~ s& , so the integrated distribution Ip(s) = J Q S Pp{s')ds' behaves as 
Ip{s) ~ s 1+/3 . A matrix U of size N yields N spacings Sj. Thus the minimal spacing s min 
occurs on average for such an argument of the integrated distribution that i/3(s m in) « 1/iV. 
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This implies that (s m in) 1+ ^ ~ ^-/N, which allows us to estimate the average minimal spacing 

(smin) ~N~rh. (3.1) 

For CUE, i.e. in the case (3 = 2 this statement is consistent with the rigorous results [TT] of 
Arous and Bourgade. As shown in Fig. 13. li the above heuristic reasoning provides the correct 
value of the exponent in dependence of the mean minimal spacing (s m i n ) on the matrix size 
N for CPE {(3 = 0), COE (J3 = 1) and CUE {fi = 2). 




Figure 3.1: Mean minimal spacing (s m ; n ) as a function of the matrix size N — 2™ for (v) 
CPE, (□) COE and (o) CUE and m = 1, 7. Symbols denote numerical results obtained for 
2 14 independent random matrices. Solid, dashed and dash-dot lines are plotted with slopes 
implied by the estimation (|3.1|) and equal to —1, —1/2 and —1/3, respectively. Linear fit to 
numerical data yields the slopes -0.98, -0.48, -0.33, respectively. 



3.2 Mean maximal spacing 

We are going to estimate the average maximal spacing (s max ) for random unitary matrices 
of the circular orthogonal ensemble. Matrix of size N yields N spacing Sj. In analogy to 
the previous reasoning we shall assume that all spacings are independent random variables 
described by the Wigncr surmise 

P(s) = \se~^'\ (3.2) 

Thus the integrated distribution I(s) = P(s')ds' reads I(s) = 1 — e^ 7 " / 4 . The maximal 
spacing s max occurs on average for such an argument of the integrated distribution function 
that 1 — I(s max ) fn l/N. This implies that e _7rSm » x / 4 w 1/N, which allows us to estimate the 
average maximal spacing, 

(Smax)cOB~- lnW - ( 3 - 3 ) 

This implies that (s max ) 2 grows with the matrix size N proportionally to ^\nN what is 
demonstrated in Fig. 13.21 
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Now let us deal with the case of the circular unitary ensemble. Here we employ the Wigner 
formula for the level spacing distribution P of a large CUE matrix P{s) = §|s 2 e~ 4s By 
the same reasoning as above we obtain an estimate that the maximal spacing s max occurs on 
average for such an argument of the integrated distribution function Z(s) = J Q S P(s')ds' that 
1 - I(w) « UN. Thus 

We change the variable u = As 2 /it and obtain w J^ 2 ^ -^u 1 ^ 2 e~ u du. Therefore, sup- 
posing s max is large we get 

N " ir"~ 

Now we take the logarithm of both sides, neglect In s max as it is of lower order than s^ ax for 
large s max , and arrive at 

(s ma x) 2 cuE ~ j lnN - (3-6) 
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■Smax e max • (3-5) 




Figure 3.2: Mean maximal spacing (s max ) as a function of the matrix size N = 2 m for a) (□) 
COE and (o) CUE, b) ( V) CPE and m = 0, 10. Symbols denote numerical results obtained 
for 2 14 independent random matrices. Solid, dashed and dash-dot lines are plotted with slopes 
implied by the estimation (3.5), (3.8), (3.9) and equal to i/n, 1 and 7r/4, respectively. Linear 
fit to numerical data yields the slopes 1.32, 0.97, 0.84, respectively. 



In the case of a Poissonian spectrum the level spacing distribution displays an exponential 
tail, P(s) ~ e~ s . Thus the integrated distribution function I(s) = J P(s')ds' behaves as 
I(s) = 1 — e _s . For a matrix of size N the maximal spacing s max occurs on average for such 
an argument that 1 — I(s max ) rj 1/N. This implies that e _Smax « 1/N and enables us to 
estimate the average maximal spacing for the circular Poisson ensemble as 

(Smax)CFB ~ IniV. (3.7) 

Analyzing estimations from equations (|3 . 3[) . p.6p and (|3.T[) gives the slopes Aqoe = — ~ 1-27, 
^■CC/b = f « 0.77 and Acpb = 1, which are comparable with numerical results Acoe ~ 1-33, 
A C ue « 0.84 and A C pb w 0.97. 
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3.3 Distribution of extremal spacings 

To study the distributions of the minimal spacing P min introduce a rescaled variable suggested 
by (3D, 

x ( IL--=A W NTT?s milL , (3.8) 

where Arp\ is a constant, in general different for CPE, COE and CUE. 

The case of the unitary ensemble was recently studied by Arous and Bourgade who 
derived the following expression for the asymptotic distribution of the minimal spacing 

-fmin(^min) = 3x m ^ n e mm , (3.9) 

in the rescaled variable x m i n = (tt /i) 2 / 3, N 1 / 3, s m \ n . This result suggests the following general 
form of the distribution of minimal spacing for all three ensembles considered labeled by the 
level repulsion parameter /3 

P^ n (x min ) := (13 + lJa^e-S. 1 , (3.10) 

which agree perfectly with the numerical data. 

The above formula has a structure F(x) := d ^j^ g~^ x \ which helps to determine the 
normalization. Numerical results show that constants read: Ami = 1 for CPE, Am = (s) for 
COE, and A {2) = (tt/3) 2/3 for CUE. 

Returning to the original variable we obtain the distributions -P^^Smin), 

P£L = A (0) Ne- N ^», (3.11) 

P£L = ZA^Ns^e-^^, (3.12) 

P^ = 3Al 2) Ns 2 min e- A U^. (3.13) 

The distributions of the minimal spacing obtained numerically for Poisson, orthogonal and 
unitary circular ensembles of random matrices of the size ./V = 100 are presented in Fig. 3.3. 
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Figure 3.3: Probability distributions P m ; n for random unitary matrices of CUE (v), COE 
(□), and CPE (o) in normal scale (fig. a)), and rescaled for b) CUE, c) COE and d) CPE. 
Symbols denote numerical results obtained for 2 17 independent matrices of size N = 100, 
while solid curves represent asymptotic predictions. 



4 Concluding remarks 

A significant and spectacular difference between the Poissonian ensemble on one side and COE 
and CUE on the other, vis. degree of "repulsion" between adjacent levels can be effectively 
analyzed in terms of distributions of the minimal spacings. We analyzed the average minimal 
spacing and postulated a general form of its distribution for the classical ensembles. For CUE 
this distribution coincides with the one recently derived by Arous and Bourgade We 
support our findings by numerical simulations showing a perfect match with the theoretical 
predictions. 
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